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Abstract 

This  report  examines  fundamental  processes  involving  the  interaction  of  nerve  agents  with  solid 
surfaces  using  accurate  electronic  structure  methods.  The  objective  is  to  describe  surface 
interactions  and  low  energy  electronic  states  accurately  and  to  identify  factors  that  affect 
desorption  energies  and  kinetics,  solvation  of  nerve  agents  by  water  and  the  spectral  features  of 
adsorbed  nerve  agents.  Understanding  exactly  how  these  molecules  interact  with  organic  films 
and  oxide  surfaces  is  the  primary  goal  of  this  work.  Initial  work  on  sarin  adsorbed  on  aliphatic 
and  aromatic  hydrocarbon  surfaces  and  on  a  calcium  oxide  surface  has  been  completed.  Studies 
of  the  solvation  of  sarin  by  water  are  underway.  The  results  should  clarify  how  water  affects 
desorption  properties  and  the  extent  of  solvation  of  the  nerve  agent  as  it  leaves  the  surface. 
Preliminary  solvation  results  are  included  in  the  report.  One  paper,  submitted  for  publication  in 
J.  Physical  Chemistry,  is  undergoing  minor  revision  to  address  question  raised  in  the  review.  A 
second  paper  on  the  solvation  studies  in  being  written. 

1.  Overview  of  Research 

In  this  work,  we  examine  fundamental  processes  involving  the  interaction  of  nerve  agents 
with  solid  surfaces  using  accurate  electronic  structure  methods.  Adsorption  on  several  types  of 
surfaces  found  in  structural  materials  are  considered  for  two  types  of  nerve  agents:  sarin,  a  nerve 
gas,  (see  Fig.  1)  and  VX  type  nerve  agents  .  Although  both  nerve  agents  are  organophosphates, 
they  differ  in  the  ester  linkage  and  in  substituent  groups.  These  differences  are  known  to  lead  to 
different  persistence  times  on  surfaces. 
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Fig.  1.  Sarin  nerve  gas.  An  organo- 
phosphate  nerve  agent  with  an  ester 
linkage  to  an  aliphatic  group. 

The  objective  of  the  work  is  to  describe  surface  interactions  and  low  energy  electronic 
states  at  high  aaccuracy  using  theory  in  order  to  understand  the  details  of  bonding  to  surfaces. 
The  plan  is  to  carry  out  the  study  so  that  factors  that  affect  desorption  energies,  desorption 
kinetics  and  the  spectral  features  of  adsorbed  nerve  agents  can  be  quantified.  Understanding 
exactly  how  these  molecules  interact  with  a  variety  of  solid  surfaces  is  the  primary  goal  of  this 
work. 

The  project  is  divided  into  three  phases: 

1)  Adsorption  of  sarin  on  hydrocarbon  and  oxide  surfaces 

2)  Solvation  of  sarin  adsorbed  on  surfaces  by  water 

3)  Study  of  the  nerve  agent  VX  and  examination  of  the  low  lying  excited  states  of  sarin 
adsorbed  on  surfaces. 

Phase  1  has  been  completed  and  Phase  2  is  nearing  completion.  The  project  benefitted 
from  the  participation  of  two  sabbatical  leave  visitors  to  the  Whitten  group  from  Athens,  Greece, 
G.  Petsalakis  and  I.  Petsalakis.  The  visitors  had  prior  experience  using  Density  Functional 
Theory  to  describe  molecule-surface  interactions  and  this  permitted  a  comparison  of  DFT  and 
our  many-electron  Cl  calculations  for  sarin  bound  to  hydrocarbon  and  oxide  surfaces. 

Appendix  I  contains  the  manuscript  describing  the  results  from  Phase  1  of  the  project. 
This  manuscript  was  submitted  to  J.  Phys.  Chem.  and  is  currently  undergoing  minor  revision  to 
clarify  questions  raised  in  the  review. 
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In  the  Phase  2  solvation  studies,  the  surfaces  are  either  initially  coated  by  water,  see  for 
example  the  snapshot  in  Fig.  2,  followed  by  introduction  of  sarin  at  a  series  of  distances  from  the 
surface,  or  the  sarin  molecule  is  solvated  first,  see  Fig.  3,  and  then  allowed  to  approach  the 
surface.  In  both  cases,  sarin  and  the  surface  compete  for  the  water  molecules  and  the  molecules 
move  optimally  to  solvate  both  sarin  and  the  surface  as  the  molecule  approaches  the  surface.  We 
follow  the  energetics  of  the  system  as  a  function  of  the  distance  of  sarin  from  the  surface 
allowing  the  water  molecules  to  move  to  their  minimum  energy  location  for  each  choice  of  sarin 
position.  The  results  should  lead  to  an  understanding  of  the  effect  of  water  solvation  on  the 
ability  of  a  surface  to  retain  sarin. 


Fig.  2.  Binding  directly  to  water  layer:  16  kcal/mol 
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Fig.  3.  Solvation  of  sarin  near  an  aliphatic  hydrocarbon  surface  by  water. 
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Appendix 

J.  Phys.  Chem.  paper 

Cl  and  DFT  Studies  of  the  Adsorption  of  the  Nerve  Agent  Sarin  on  Surfaces 

Brian  N.  Papas,3  loannis  D.  Petsalakis,ab  Giannoula  Theodorakopoulosa,b  and  Jerry  L. 

Whitten3 

“Department  of  Chemistry,  North  Carolina  State  University,  Raleigh,  North  Carolina,  27695, 
USA. 

bTheoretical  and  Physical  Chemistry  Institute,  The  National  Hellenic  Research  Foundation,  48 
Vass.  Constantinou  Ave.  Athens  116  35  Greece. 

ABSTRACT:  A  theoretical  study  is  presented  on  the  adsorption  of  the  sarin  molecule,  a 
nerve  agent,  on  three  model  solid  surfaces:  aliphatic  (graphane),  aromatic  (graphene),  and 
ionic  (CaO).  Calculations  are  carried  out  using  accurate  electronic  structure  methods  based 
on  Configuration  Interaction  (Cl)  as  well  as  complementary  Density  Functional  Theory 
(DFT)  and  Time  Dependent  DFT  calculations.  The  objective  is  to  describe  surface 
interactions  accurately  in  order  to  identify  factors  that  affect  adsorption  of  sarin.  Potential 
energy  curves  are  calculated  and  compared  between  surface  types.  Computed  Cl  binding 
energies  to  graphane,  graphene,  and  calcium  oxide  surfaces  are  2.4,  5.2,  and  13.2  kcal/mol, 
respectively.  The  corresponding  DFT  binding  energies  are  6.4,  4.8,  and  18.8  kcal/mol, 
respectively.  Excited  states  of  free  sarin  as  well  as  sarin  bound  to  the  surfaces  are  examined 
using  TDDFT  and  Cl.  Low-lying  excited  states  of  the  adsorption  complexes  involving 
excitations  from  the  surface  to  sarin  are  calculated. 

1.  INTRODUCTION 

In  this  work,  we  examine  fundamental  processes  involving  the  interaction  of  sarin 
(isopropyl  methylphosphonoflouridate),  a  nerve  agent  molecule  depicted  in  Fig.  1,  with  aliphatic, 
aromatic,  and  polar  surfaces  using  accurate  electronic-structure  methods  of  configuration 
interaction  (Cl)  and  density  functional  theory  (DFT).  Sarin  is  a  member  of  a  class  of 
organophosphate  (OP)  nerve  agents  that  differ  in  ester  linkage  and  substituent  groups.  The 
objective  of  this  work  is  to  describe  surface  interactions  at  high  accuracy  using  theory  in  order  to 
understand  the  details  of  bonding  of  sarin  to  common  surfaces.  Detailed  insight  into  the  surface 
chemistry  of  chemical  warfare  agents  (CWAs)  is  critical  to  the  rational  design  of  advanced  filters 
and  decontamination  strategies  for  these  extremely  toxic  compounds. 
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Figure  1.  Sarin  nerve  gas.  An  organo-phosphate  nerve  agent  with  an  ester  linkage  to  an 
aliphatic  group.  Coloring  scheme  as  follows:  C  (gray),  H  (white),  O  (red),  P  (orange),  F  (blue). 

Most  nerve  agents  act  by  disrupting  the  enzymatic  degradation  of  the  neurotransmitter 
acetylcholine.  The  organophosphates  are  known  to  bind  nearly  irreversibly  to 
acetylcholinesterase,  thereby  blocking  the  channel  for  acetylcholine  decomposition  and 
recycling.  The  resulting  buildup  of  acetylcholine  over-stimulates  nerve  action  leading  to  fatal 
contractions  of  muscles. 

The  same  molecular  characteristics  that  cause  tight  binding  to  acetylcholinesterase  may 
be  involved  in  bonding  to  structural  materials.  It  is  useful  to  consider  first  the  interactions 
qualitatively.  Strong  interactions  involving  the  phosphate  region  and  ester  linkages  can  induce 
ionic  or  polar  interactions  with  a  solid  surface,  and  presumably  this  would  leave  the  molecule 
free  to  later  desorb  as  an  intact  molecule.  Strong  covalent  bonding,  however,  may  disrupt  the 
molecular  structure  and  leave  the  agent  inactive.  Substituent  groups  can  give  rise  to  a  stabilizing 
van  der  Waals  bonding  with  any  surface,  and  if  only  this  type  of  bonding  were  involved,  one 
would  expect  that  the  adsorbed  molecule  could  be  easily  removed  from  the  surface. 
Understanding  quantitatively,  at  the  molecular  level,  the  interactions  of  sarin  with  the  different 
surfaces  and  the  excited  states  of  the  adsorption  complexes  are  the  main  objectives  of  the  present 
investigation. 
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Figure  2.  Model  surfaces  employed  in  the  Cl  calculations.  Left-to-right:  graphane,  graphene, 
calcium  oxide.  Covalent  radii  used  for  atom  size.  Coloring  scheme  as  follows:  C  (gray),  H 
(white),  O  (red),  Ca  (yellow). 

Two  types  of  hydrocarbon  surfaces  and  calcium  oxide  are  adopted  as  models: 

1.  Graphane  -  a  surface  with  CH  dipoles  and  low  polarizability. 

2.  Graphene  -  polycyclic  aromatics  show  intermediate  reactivity  and  the  potential 
for  large  dipole  and  charge  induced  polarizations. 

3.  Calcium  oxide  -  the  ionic  solid  limit  where  strong  interactions  with  the  phosphate 
and  ester  regions  of  the  nerve  agents  are  anticipated.  Calcium  oxide  surfaces  have 
similar  ionic  charge  distributions  to  calcium  carbonates.  Carbonates  and  silicates 
are  major  components  of  structural  materials. 

A  great  deal  of  experimental  and  computational  work  has  been  devoted  to  adsorption  and 
interaction  of  OP  warfare  agents  or  stimulants  with  different  surfaces,  including  silica,1’ 2 
carbonaceous  nanoporous  materials,  and  metal-oxide  surfaces.3"8  To  the  author’s  knowledge, 
adsorption  of  sarin  on  the  surfaces  of  the  present  work,  graphane,  graphene  and  CaO,  has  not 
been  previously  considered,  while  there  exists  a  previous  theoretical  investigation  of  sarin 
physisorbed  as  well  as  chemisorbed  on  model  MgO  surfaces.6  The  structure  of  gas  phase  sarin 
has  also  been  investigated  extensively9"17  by  spectroscopic  methods  and  by  electronic-structure 
calculations.  Furthermore,  dissociation  of  sarin  in  the  lowest  triplet  and  singlet  excited  states  has 
been  investigated16  in  response  to  a  report  of  destructive  elimination  of  sarin  with  the  aid  of  UV 
irradiation.14’ 15 

2.  THEORETICAL  METHODS 
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Two  types  of  electronic  structure  approaches  are  employed:  the  first  one  involves 
configuration  interaction  along  with  a  simplex  optimization  method  for  the  geometry 
optimizations.  The  second  approach  involves  density  functional  theory  (DFT)  calculations  for 
the  ground  electronic  state  and  time  dependent  DFT  (TDDFT)  calculations  for  excited  electronic 
states,  with  details  given  below. 

2.1.  Many-electron  Cl  Theory 

The  adsorbate-surface  systems  will  be  described  by  a  configuration  interaction 
embedding  method  that  is  designed  to  give  an  accurate  many-electron  description  of  the 
adsorbate-surface  region  for  both  ground  and  excited  electronic  states.  Many  applications  to 
adsorbates  on  metals  and  oxides  have  been  reported  previously  and  details  can  be  found  in  Ref. 
18-21.  In  the  case  of  small  particles,  there  is  no  need  to  make  approximations  in  the  coupling  of 
the  local  region  to  the  bulk  and  in  this  case  “embedding”  refers  to  the  creation  of  an  electronic 
subspace  that  is  treated  by  configuration  interaction.  A  brief  summary  of  the  theory  is  given 
below. 

Calculations  are  carried  out  for  the  full  electrostatic  Hamiltonian  of  the  system 


N  Q  y  N 

h  =  Z  hv;-+  £(—)]+£ 

i  k  rik  i<j 
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Wavefunctions  are  constructed  by  self-consistent-field  (SCF)  and  multi-reference 
configuration  interaction  (Cl)  expansions, 

'V  =  (n!  f  det  (Zl‘  Z‘  ■  ■  ■  X*  )  =  2>t<t>t 

k  k 

Single  and  double  excitations  from  an  initial  representation  of  the  state  of  interest,  Or,  are 
carried  out  to  create  a  small  Cl  expansion, 

TV  =  Or  +  /„jjki  r ij^ki  Or  =  Em  Om 

Configurations,  Ok,  are  retained  if  the  interaction  with  Or  satisfies  a  relatively  large  threshold 
condition 


l<  Ok  I  H  I  Or  >1  2  /  IEk  -  Er  I  >  10  4  a.u. 
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Next,  we  refine  the  description  of  the  state  by  generating  a  large  Cl  expansion,^,  by  single  and 
double  excitations  from  all  important  members  of  TV  (those  with  coefficient  >  0.05)  to  obtain 

*Pr  =  VF/  +  ^mPik  A-iktn  r^kOm  +  2)jkl  ^ijklm  Tjj_>kl  CbmJ 

where  <J>m  is  an  important  member  of  TV-  The  additional  configurations  are  generated  by 
identifying  and  retaining  all  configurations,  Ok,  that  interact  with  TV  such  that 

l<  Ok  I  H  I  TV  >1  2  /  IEk  -  Er  I  >  3x10  7  a.u. 

For  this  small  threshold,  typically  4xl04  configurations  occur  in  the  final  Cl  expansion,  and  the 
expansion  can  contain  single  through  quadruple  excitations  from  an  initial  representation  of  the 
state  Or.  Contributions  of  configurations  not  explicitly  retained  are  estimated  using  perturbation 
theory. 

In  order  to  accelerate  the  convergence  of  the  Cl  expansion  results  with  respect  to  active 
space  size  for  the  excited  state  calculations,  the  molecular  orbitals  of  the  lowest  triplet  state  for 
each  system  are  first  put  through  a  unitary  transformation.  This  transformation  creates  a  set  of 
occupied  and  virtual  orbitals  which  have  maximal  exchange  interaction  with  the  occupied 
orbitals  within  the  active  space.  The  primary  function  of  this  is  to  allow  accurate  description  of 
the  states  of  interest  in  the  presence  of  less  accurately  described  lower-lying  states.  The  basis  of 
this  argument  lies  in  the  way  other  configurations  are  generated  from  the  states  of  interest.  For  a 
given  excited  state,  TV  allowing  all  single  and  double  excitations  from  its  important 
configurations  also  generates  all  lower  lying  states  that  have  significant  interaction  with  TV 
Since  these  configurations  are  in  principle  included  in  the  diagonalization,  it  follows  that  the 
excited  state,  TV  is  orthogonal  to  and  non-interacting  with  lower  energy  states.  Thus  the 
variational  theorem  for  excited  states  is  satisfied. 

For  the  Cl  computations,  the  interaction  of  sarin  is  described  using  fixed-geometry 
surfaces  sufficiently  large  enough  to  ensure  that  sarin  would  remain  fully  over  the  surface  and 
not  interact  appreciably  with  boundary  atoms.  These  surface  models,  depicted  in  a  subsequent 
figure,  correspond  to  C22H12  for  graphene,  C22H34  for  graphane,  and  Ca^Ois  for  calcium  oxide. 

Compact,  but  still  flexible  bases  were  used  in  this  study,  and  can  roughly  be  characterized 
as  double-zeta  plus  polarization.  These  are  described  as  (10s5pld)— >{4s2pld]  for  C  and  O, 
(12s8pld)— >[6s4pld]  for  P,  (10s5p)— >[4s2p]  for  F,  and  (5slp)— >[2slp]  for  H.  For  the  graphane 
surface  only  the  pz  orbital  was  added  to  each  H  atom  to  ensure  sufficient  bond  polarization. 
Graphene  hydrogen  atoms  did  not  use  any  p  functions,  while  the  carbon  atoms  had  additional 
diffuse  s  and  pz  functions  added  with  exponents  of  0.08.  For  calcium  oxide  the  surface  should 
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show  little  reaction  to  the  addition  of  sarin  and  so  a  (22s  12p)— »[5s2p]  basis  for  Ca  and  a 
(1 2s5p) — >[2s  1  p]  basis  for  O  were  used.  This  basis  was  generated  by  minimizing  the  total  energy 
of  Ca4C>4,  and  is  given  in  the  supplemental  information. 

2.2.  Simplex  Optimization  Method 

The  Cl  geometry  optimizations  are  carried  out  by  the  Nelder-Mead  simplex  procedure.23 
The  optimization  allows  variation  of  the  geometry  of  adsorbates  on  the  surfaces  of  interest.  In 
the  simplex  procedure,  the  parameters  that  define  the  geometry  are  taken  as  components  of  a  so- 
called  vertex  or  n-tuple.  The  energy  is  calculated  for  that  vertex  and  a  systematic  variation  of  the 
parameters  is  carried  out  in  a  way  that  eliminates  the  worst  vertex  after  each  iteration.  The 
procedure  is  widely  used  in  applied  mathematics  and  is  surprisingly  robust.  It  has  an  advantage 
in  that  energy  gradients  are  not  required  and  the  energy  minimization  can  be  carried  out  directly 
at  the  Cl  level  of  theory  since  only  a  computer  code  capable  of  computing  the  energy  is  needed. 
The  latter  point  is  important  when  electron  correlation  has  a  significant  effect  on  the  geometry, 
e.g.,  stretched  bonds  or  pre-dissociation.  Geometry  optimizations  at  the  Cl  level  allowed  all 
degrees  of  freedom  of  sarin  to  vary,  while  keeping  the  model  surfaces  locked  at  a  standard 
geometry. 

2.3.  DFT  and  TDDFT  Calculations 

Complementary  to  the  Cl  calculations  described  in  Section  2.1  DFT24  and  TDDFT25 
calculations  employing  the  M062X  functional  ’  and  two  types  of  basis  sets,  6-3  lG(d,p)  and  6- 
31  lG+(d,p)  (i.e.  including  diffuse  functions)  have  been  carried  out  on  sarin  adsorbed  on 
graphene,  graphane  and  CaO,  as  in  the  Cl  calculations  but  using  slightly  different  models  for  the 
surfaces.  The  smaller  basis  set  was  employed  for  larger  models  of  the  surfaces,  cf.,  structures  (a) 
and  (d)  for  graphene  (C66H22)  and  graphane  (C66H88),  respectively,  and  structure  (g)  for  CaO  in 
Figure  3.  The  larger  basis  set  was  employed  for  smaller  models  of  the  graphene  (C28H14)  and 
graphane  (C28H42)  surfaces,  (b)  and  (e),  respectively,  in  Figure  3,  and  also  for  the  CaO  (Ca^Ois) 
surface,  (h).  Binding  energies  of  sarin  on  the  three  surfaces  have  been  determined,  taking  into 
account  the  basis  superposition  error.  In  addition,  excited  states  and  absorption  spectra  of  the 
adsorption  complexes  of  sarin  on  graphene  and  graphane  have  been  calculated  and  compared  to 
the  absorption  spectrum  of  free  sarin  obtained  by  TDDFT/M062X/6-31 1+G(d,p).  For  further 
comparison,  the  absorption  spectrum  of  isolated  sarin  has  been  calculated  by  TDDFT  employing 
nine  additional  functionals  (specified  in  Figure  4)  and  the  6-311+G(d,p)  basis  set.  Corresponding 
data  has  been  generated  by  Cl  calculations  for  free  sarin  as  well  as  the  adsorption  complexes. 

Geometry  optimization  of  the  adsorption  complexes  involved  first  the  optimization  of  the 
surface  and  subsequently  optimization  of  the  adsorbed  sarin.  In  the  case  of  CaO,  inclusion  of  5 
atoms  on  the  surface  in  the  second  step  of  the  optimization  resulted  in  2.3  kcal/mol  increased 
binding  energy. 
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Finally,  DFT/M062X  calculations  have  also  been  carried  out  for  sarin  on  MgO  (Mg^Ois) 
(cf.  (h)  in  Figure  3),  for  comparison  with  previous  work  on  this  system.6  For  these  calculations 
basis  sets  6-31G(d,p)  and  6-311G(d,p)  have  been  employed.  Addition  of  diffuse  functions  was 
not  possible  as  the  calculations  did  not  converge.  For  comparison,  the  previous  work  employed 
6-31G(d)  basis  set  with  DFT/B3LYP  and  MP2  calculations  on  different  models  of  sarin  on 

i. : 

MgO.  All  the  DFT  and  TDDFT  calculations  of  the  present  work  have  been  carried  using 
Gaussian  09.28 
3.  RESULTS 

3.1.  Sarin  Geomerties 

Four  conformers  of  sarin  were  optimized  at  the  Cl  level.  The  minimum  energy  structure 
of  sarin  used  in  the  present  study  is  depicted  in  Figure  1 .  This  is  the  geometry  labeled  sarin  II 
reported  in  Ref.  17.  The  energy  of  this  geometry  is  less  than  1  kcal/mol  above  the  minimum 
energy  geometry  reported  in  Ref.  17,  and  thus  for  the  purposes  of  this  research  is  considered 
unimportant.  The  remaining  conformers  differed  by  rotation  of  the  aliphatic  and  phosphate 
regions  and  changes  in  the  central  C-O-P  bond  angle. 

3.2.  Binding  Geometries 
3.2.1.  Overview 

In  Figure  3  the  structures  of  the  adsorption  complexes  obtained  by  DFT  and  Cl  are 
given;  in  Table  1  the  calculated  binding  energies  are  collected.  These  results  will  be  discussed 
below  starting  with  the  comparatively  simplest  case  of  sarin  on  CaO. 
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(a)  Graphene,  DFT 


a=3.13  A 
b=3.11  A 


(b)  Graphene,  DFT  (c)  Graphene,  Cl 


A 

A 

A 

A 


(d)  Graphane,  DFT 


(e)  Graphane,  DFT 


(f)  Graphane,  Cl 


(g)  CaO,  DFT 


(h)  MgO,  DFT 


(i)  CaO,  Cl 


Figure  3.  Optimum  geometries  of  the  adsorption  complexes  of  sarin  on  graphene,  (a)  and  (b) 
obtained  by  DFT  and  (c)  by  Cl,  on  graphane,  (d)  and  (e)  by  DFT  and  (f)  by  Cl,  on  CaO,  (g)  by 
DFT  and  (i)  by  Cl  and  on  MgO  (h)  by  DFT,  cf.  text. 

3.2.2.  Sarin  on  CaO 

Figures  3g — i  show  optimal  binding  geometries  of  sarin  to  oxide  surfaces.  The  binding 
appears  to  be  a  straight-forward  electrostatic  interaction.  The  P  atom  of  sarin  aligns  with  an  O 
atom  of  the  surface,  and  the  F  and  terminal  O  atoms  align  over  Ca  atoms,  creating  pairs  of 
opposite  charges.  The  center  of  mass  of  sarin  is  4.04  A  above  the  surface  in  the  Cl  calculations. 
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In  all  calculations,  the  Ca-O/F  distance  is  on  the  order  of  2.4  A.  All  calculations  are  in  good 
agreement  on  the  binding  energy  (Table  1),  with  all  results  in  the  range  of  10-20  kcal/mol. 
Relaxation  of  the  CaO  lattice  during  binding  was  found  to  have  negligible  impact  on  the  binding 
energy  using  DFT.  For  comparison,  reported  values  of  binding  energy  of  sarin  on  MgO, 
calculated  by  B3LYP/6-31G(d)  and  a  large  model  (Mgi60i6)  is  2.9  kcal/mol  while  that  obtained 
by  MP2/6-31G(d)  and  a  small  cluster  (Mg404)  is  50  kcal/mol.6  It  should  be  noted  that  in  the 
present  work  the  MgO  was  modelled  with  M062X  on  a  single  layer  of  Mgi50i5,  not  two  layers 
as  in  the  previous  work.  As  shown  in  Table  1,  significant  binding  of  20.0  kcal/mol  is  obtained 
for  sarin  on  MgO,  consistent  with  the  expectation  of  stronger  binding  of  sarin  to  MgO  than  CaO. 

Table  1.  Calculated  binding  energy  (BE)  of  sarin  at  different  surfaces  by  DFT  and  CL 


Surface  (structure  of  Figure  3)  /basis  set 

BE  (kcal/mol) 

BSSE-corrected 

BE  (kcal/mol) 

Graphene  (a)/  6-31G(d,p) 

6.7 

- 

Graphene  (b)/  6-31G(d,p) 

6.4 

2.6 

Graphene  (b)/  6-311+G(d,p) 

8.2 

4.8,  5.2* 

Graphane  (d)/  6-31G(d,p) 

9.0 

- 

Graphane  (e)/  6-31G(d,p) 

7.9 

4.9 

Graphane  (e)/  6-311+G(d,p) 

8.5 

6.4,  2.4* 

CaO(g)/  6-31G(d,p) 

33.9 

16.1 

CaO(g)/  6-311+G(d,p) 

22.3 

18.8,  13.2* 

MgO  (h)/6-31G(d,p) 

32.8 

19.3 

MgO  (h)/6-311G(d,p) 

31.8 

20.0 

^Corresponding  Cl  value 


3.2.3.  Sarin  on  graphene 

Unlike  the  results  for  CaO,  DFT  and  Cl  give  different  binding  geometries  for  sarin  on 
graphene  (Figure  3a-c).  Specifically,  DFT  models  predict  binding  through  the  sarin  oxygen  and 
fluorine  atoms,  while  Cl  predicts  a  flipped  orientation,  cf.  Figure  3b  vs  3c.  Consistency  between 
the  M062X/6-31G(d,p)  results  for  large  (Fig.  3a)  and  middle-sized  (Fig  3b)  graphene  surfaces 
suggests  that  the  M062X/6-31 1+G(d,p)  binding  energy  of  4.8  kcal/mol  should  be  taken  as 
accurate.  Despite  the  difference  in  orientation,  there  is  good  agreement  with  the  Cl  computed 
binding  energy  of  5.2  kcal/mol.  Cross-examinations  of  the  two  orientations  with  the  two  levels 
of  theory  indicate  that  the  energy  difference  between  the  two  orientations  is  small  (about  0. 1-0.4 
kcal/mol,  depending  on  the  model),  with  the  optimum  geometry  being  dependent  on  the 
particulars  of  the  method  used.  It  should  be  noted  that  in  both  orientations  the  dipole  moment  of 
sarin  is  directed  perpendicular  to  the  graphene  surface. 
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3.2.4.  Sarin  on  graphane 

Of  the  three  surfaces  modelled,  graphane  exhibits  the  greatest  discrepancy  between  the 

O 

DFT  and  Cl  results.  From  the  DFT  results,  sarin  binds  with  its  oxygen  atom  down  and  2.3  A 
from  a  nearest  surface  hydrogen  atom,  and  its  fluorine  atom  pointed  away  from  the  surface  (Fig. 
3d,e).  Predicted  binding  is  on  the  order  of  6  kcal/mol  -  slightly  greater  than  that  of  sarin  on 
graphene.  The  Cl  computations  suggest  a  binding  with  the  sarin  oxygen  and  fluorine  atoms  at 
the  surface  (Fig.  3f),  but  a  binding  of  only  2.4  kcal/mol.  Investigation  of  the  orientation  of  sarin 
on  graphane  obtained  by  Cl  employing  the  Cl  model  (C22H34),  as  well  as  the  larger  (C28H42) 
model,  using  the  DFT/M062X  method  show  that  for  the  smaller  Cl  model  only  the  Cl  structure 
(cf.  (f)  of  Figure  3)  is  obtained  whereas  with  the  larger  model  the  DFT  structure,  cf.  (e)  in  Figure 
3,  is  favored.  Therefore  there  is  a  significant  effect  of  the  size  of  the  model  employed  on  the 
calculated  orientation  of  the  adsorption  complex. 

3.3.  Excited  States 

A  great  deal  of  experimental  and  theoretical  work  has  been  devoted  to  the  determination 
of  the  lowest  energy  structure  of  sarin  and  other  warfare  agents  in  the  ground  electronic  state,  but 
the  excited  states  have  not  been  studied  very  extensively.5  For  sarin,  a  low  intensity  absorption 
is  reported  with  the  threshold  electronic  excitation  occurring  in  the  vacuum-UV  (~7.1  eV  or 
higher)  with  varying  theoretical  estimates,  eg.  7.38  eV  obtained  by  TDDFT/B3LYP/6- 
31 1+G(2df,p)  and  9.9  eV  by  MP2/6-31+G(d)  calculations.9  Here  too  we  find  a  varying 
excitation  energy  for  the  lowest  singlet  excited  state  of  free  sarin,  depending  on  the  functional 
employed  in  TDDFT/6-31 1+G(d,p)  calculations  .  In  Figure  4  the  absorption  spectra  of  free  sarin 
obtained  by  TDDFT  calculations  and  ten  different  functionals,  offered  in  Gaussian  09,  are 
shown.  The  lowest  onset  for  the  electronic  excitations  is  obtained  with  the  PBEPBE  functional 
at  6.6  eV  and  the  highest  with  M062X  at  8.1  eV.  Generally  low  oscillator  strengths  are 
calculated  with  all  functionals  employed.  Cl  calculations  on  the  excited  state  employing  the 
basis  set  previously  described  lead  to  an  excitation  energy  of  9.1  eV,  i.e.,  comparable  to  an 
earlier  MP2/6-31+G(d)  result.9  However,  inclusion  of  a  single  diffuse  s  function  (exponent  of 
0.01)  in  the  phosphorus  atom  basis  set  allows  an  accurate  Cl  estimate  of  the  lowest  Rydberg 
excitations.  These  come  in  at  8.0  eV  and  8.5  eV,  corresponding  to  excitations  from  the 
oxygen/fluorine  lone-pair  manifold  into  a  Rydberg  state.  This  is  in  good  agreement  with  the 
present  CAM-B3LYP  and  M062X  TDDFT  results  on  the  spectrum  of  free  sarin,  cf.  Figure  4 
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Figure  4.  Adsorption  spectra  of  free  sarin  calculated  by  TDDFT  employing  different 


functionals,  as  indicated. 


Next,  the  absorption  spectra  of  sarin  adsorbed  on  graphene  and  graphane  were  considered 
using  TDDFT/M062X  and  the  6-31 1+G(d,p)  basis  set.  For  the  graphene-sarin  system  the  lowest 
50  roots  were  calculated  in  the  region  1.98  eV  -  6.19  eV,  none  of  which  correspond  to  the  sarin- 
sarin  excitation.  The  lowest  transitions  correspond  to  graphene-graphene  excitations.  There  also 
exist  a  number  of  charge-transfer  type  excitations,  from  occupied  orbitals  of  graphene  to 
unoccupied  orbitals  of  sarin.  For  example,  states  calculated  at  3.37  and  4.2  eV  involve 
excitations  into  a  Rydberg  unoccupied  orbital  of  sarin.  Other  such  states  are  found  at  5.5,  6.08, 
6.18  and  6.19  eV.  The  calculations  on  graphane-sarin  involved  the  lowest  20  roots.  In  this  case 
the  lowest  unoccupied  orbital  or  LUMO  is  a  Rydberg  orbital  of  sarin  and  is  involved  in  several 
excited  states  of  the  complex.  In  particular,  the  lowest  excited  state  at  6.4  eV,  which  is 
characterized  by  the  HOMO— >LUMO  excitation,  corresponds  to  an  excitation  from  graphane  to 
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sarin.  As  was  the  case  for  graphene-sarin,  the  calculations  did  not  reach  any  sarin-sarin  excited 
states.  It  should  be  noted  that  the  above  calculations  on  the  excited  states  of  the  adsorption 
complexes,  just  as  in  the  case  of  free  sarin,  required  inclusion  of  diffuse  functions  in  the  basis 
set.  With  the  6-31G(d,p)  basis  set  the  first  25  roots  of  the  CaO-sarin  calculation  correspond  to 
excitations  involving  only  the  surface.  Calculations  with  the  6-31 1+G(d,p)  basis  on  the  CaO- 
sarin  complex  were  not  practical. 

Cl  Computations  in  which  the  dominant  component  of  the  excitation  was  restricted 
to  be  into  the  Rydberg  orbital  of  sarin  were  used  to  determine  the  transition  energy  of  the  lowest 
lying  surface-to-sarin  electron  transfer  excitation.  For  graphane,  the  lowest  is  at  7.0  eV;  the 
remainder  are  above  the  free  sarin  excitation.  The  numerous  low-lying  virtual  orbitals  on 
graphene  proved  a  significant  challenge  for  isolating  the  charge-transfer  excitation,  and  the  best 
estimate  places  it  near  5.3  eV.  Excitations  from  the  central  portion  of  the  CaO  model  surface 
(excitations  from  orbitals  localized  on  the  edges  of  the  model  were  disallowed)  provided  a  set  of 
ten  electron-transfer  possibilities  ranging  in  transition  energy  from  6.2  eV  to  8.5  eV.  Therefore, 
with  both  methods,  the  existence  of  charge-transfer  excited  states  at  energies  lower  than  the 
sarin-sarin  excitation  is  indicated.  Similar  observations  have  been  reported  for  sarin  and  other 
organophosphorus  warfare  agents  adsorbed  on  a  Y-AI2O3  surface.5 
4.  CONCLUSIONS 

Our  results  indicate  that  the  polarized  phosphate  region  of  the  sarin  molecule  plays  a 
critical  role  in  its  binding  to  aliphatic,  aromatic,  and  oxide  surfaces.  Specifically,  the  Lewis-base 
portion  acts  as  the  binding  agent.  For  the  surfaces  considered  in  the  present  study,  only  those 
that  are  ionic  (CaO)  or  have  large  polarizability  (graphene)  give  strong  binding  for  sarin. 

Binding  to  aliphatic  surfaces  is  weak.  Figure  5  depicts  Cl  potential  curves  for  interaction  of 
sarin  with  the  three  model  surfaces. 

E/kcalmol-1 


Figure  5.  Potential  curves  for  sarin  interaction  with  three  model  surfaces.  Infinite 
separation  limit  is  free  sarin  in  minumum  energy  geometry  and  bare  surface.  The  z-axis 
corresponds  to  the  height  of  the  sarin  center-of-ma 
ss  above  the  highest  atom  of  surface. 
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It  is  concluded  that 

•  Sarin  binds  to  graphane,  graphene,  and  calcium  oxide  surfaces  by  2.4,  5.2,  and  13.2 
kcal/mol,  respectively  according  to  Cl,  with  corresponding  DFT  binding  energies  6.4, 
4.8,  and  18.8  kcal/mol,  respectively. 

•  Free  sarin  has  an  excitation  into  a  Rydberg  state  at  about  8  eV. 

•  Each  of  the  three  surfaces  studied  is  capable  of  transferring  an  electron  to  sarin  at  or 
below  the  sarin  Rydberg  excitation.  Both  DFT  and  many-electron  computational 
methods  gave  similar  results,  with  TDDFT  capable  of  determining  a  greater  total  number 
of  roots,  and  Cl  capable  to  giving  accurate  results  for  specific  states  of  interest. 
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1.  Overview  of  Tasks 

In  the  current  funded  period,  the  oxygen  reduction  reaction  (ORR)  in  the  cathode  of  solid  oxide 
fuel  cells  (SOFC)  by  molten  carbonate  (MC)  was  further  explored  using  DFT  methods.  Previous 
studies  indicate  that  the  formation  of  peroxocarbonate  (CO42  )  is  a  key  process  in  ORR  and  the 
formed  reaction  intermediate  serves  as  a  carrier  of  oxygen  in  MC.  Simply,  CO42'  can  be  treated 
as  an  atomic  oxygen  attached  to  one  carbonate.  The  migration  of  oxygen  in  MC  has  been  studied 
by  DFT,  showing  the  oxygen  in  the  form  of  CO42'  has  lower  reduction  potential  than  oxygen 
molecule  and  can  migrate  in  MC  with  low  barrier  of  energy. 

n  +co2~ _ >co2~ 

w2 r^iv3  7  ^^5  (MC) 

CO2-  +  CO2  - >209;- 

In  this  part,  we  calculated  the  formation  of  CO42'  as  a  two  step  process  (Equations  above).  First, 
molecular  oxygen  captured  by  a  carbonate  and  forms  CO52",  then  another  carbonate  will  be 
reacted  with  CO52'  to  form  two  CO42'.  Transition  states  were  located  on  the  potential  energy 
surface  (PES)  and  intrinsic  reaction  coordinate  (IRC)  calculations  show  the  reaction  process  as 
indicated  by  the  equations  above. 

The  DFT  calculated  results  will  be  applied  to  the  studies  by  the  configuration  interaction  (Cl) 
method.  In  particular,  the  optimized  structures  will  be  used  for  Cl  calculations.  The  Cl 
wavefunction  will  then  be  used  to  analyze  the  chemical  process  and  compared  to  the  DFT 
results.  The  direct  comparison  will  validate  the  accuracy  of  DFT  and  assure  the  use  of  right 
hybrid  functional(s)  when  DFT  is  applied  to  such  systems.  More  importantly,  the  Cl  method 
using  charge  localization  method  will  allow  to  examine  the  bonding  and  energy  of  charged 
molecules  in  MC  or  on  cathode  surfaces. 

2.  Results  and  Discussions 
2.1  Formation  of  CO42'  in  MC 

The  binding  of  oxygen  to  carbonate  forms  CO52',  the  optimized  structure  of  CO52'  (Figure  1)  is 
very  close  to  that  in  gas  phase  reported  and  only  negligible  change  was  detected.  This  CO52"  will 
then  react  with  another  carbonate  to  form  CO42'  as  shown  in  the  equations  below. 

For  step  two,  a  transition  state  (Figure  1)  was  located  on  the  potential  energy  surface  (PES)  for 
each  type  of  MC.  At  TS,  the  0-0  distance  is  increased  to  about  2.0A,  indicative  of  a  breaking  O- 
O  bond.  IRC  calculations  (Figure  2)  show  the  energy  barriers  and  confirm  the  reaction  path. 
More  details  will  be  included  in  the  manuscript  under  preparation.  The  shift  of  one  oxygen  to  the 
other  carbonate  causes  the  dissociation  of  oxygen  attached  to  one  carbonate,  which  makes  two 
independent  CO42',  which  then  serves  as  oxygen  carrier. 
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(A)  02-(Li2C03)4 


(B)  02-(Na2C03)4 


TS 


Figure  1:  Optimized  structures  of  CO52',  TS,  and  2CO42'  in  MC 
(0-0  distance:  -1.4 A  in  COs2 ,  -2.0A  in  TS,  and  >3.0A  in  2C042) 


21 


2-  2-  2- 

Figure  2:  PES  of  C05  +  C03  2CC>4  (Top:  Li2C03,  Middle:  Na2C03,  Bottom:  K2C03) 
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The  energy  barrier  of  oxygen  dissociation  to  form  CO42'  is  estimated  to  be  49,  28,  41  kcal/mol  in 
Li,  Na,  K  MC,  respectively.  This  finding  needs  confirmation  form  experiments  and  also  further 
modeling  with  different  methods  will  be  performed  to  confirm  if  sodium  carbonate  will  benefit 
the  ORR  process  more  than  lithium  and  potassium  carbonate  salts. 

2.2  New  Oxygen  Reduction  Mechanism  by  DFT  Modeling 

Combining  the  experimental  and  computational  results,  the  evidences  point  to  that  oxygen 
reduction  is  through  the  formation  of  CO42'  in  MC.  This  active  oxygen  form  will  accept  electrons 
and  release  an  oxide  and  carbonate.  Previous  studies  show  the  oxygen  can  easily  shift  between 
two  neighboring  carbonates,  which  means  the  oxygen  diffusion  in  MC  is  very  effective.  The  full 
reaction  path  is  displayed  in  Figure  3. 


Lfthfum:  24.3  Kcal/mol 
Sodium:  3.1  Kcal/mol 
Potassium:  16.4  Kcal/|wdf 


^  Transition  State 


Oz  +  2C032" 


Pseudo- 

Activation 

energy 


Lithium:  23.3  Kcal/mol 

abdlufn:  IS. 2  Kcal/mol 

Potassium:  25.0  Kcal/mol 


2C042 


Figure  3:  Formation  of  CO42'  by  Oxygen  in  MC 


(The  energy  barrier  for  overall  reaction  is  calculated  to  be 
24.3  kcal/mol,  3.1  kcal/mol  and  16.4  kcal/mol) 
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2.3  Raman  and  DFT  study  of  pyrocarbonate 

2.3.1  In-situ  Raman  spectroscopic  investigation 

The  Li2C03  and  Na2C03  eutectic  mixture  (52:48  in  mol%)  was  first  synthesized  in  air  at  650°C. 
After  a  2-h  hold,  the  melt  was  then  quenched  to  room  temperature,  followed  by  breaking  it  into 
fine  particles  by  ball  milling.  Thus  prepared  powders  were  then  packed  into  a  gold  crucible  that 
was  subsequently  loaded  into  a  high  temperature  stage  (Linkam  TS1500,  0°C~1500°C).  The 
temperature  and  gas  were  controlled  by  a  system  controller  (Linkam  PE95).  The  Raman  spectra 
were  recorded  with  a  LabRam/HR  confocal  Raman  system  (LabRam  Invers,  Horiba  Jobin-Yvon) 
with  a  He-Ne  laser  operated  at  632.8  nm.  Since  the  position  of  the  thermal  couples  in  the  high 
temperature  stage  is  located  outside  the  crucible,  the  actual  and  the  controlled  temperatures  of 
the  MC  are  different.  The  melting  point  of  the  of  (Li/Na)2CC>3  at  490°C  was  used  to  calibrate  the 
actual  temperature.  The  scattering  Raman  spectra  were  collected  in-situ  from  the  carbonate  as  a 
function  of  temperature  (in  the  range  of  RT-  600°C)  and  atmosphere  (in  N2,  air,  C02). 

2.3.2  DFT  modeling 

DFT  calculations  were  performed  at  the  B3LYP/6-31G(d)  level  using  Gaussian09  suite  of 
quantum  programs.  The  geometry  of  Li2C2Os  and  Na2C2Os  were  optimized  first,  and  the 
vibrational  frequencies  were  then  obtained  from  analytic  second  derivatives  using  a  harmonic 
oscillator  model.  In  addition,  Raman  intensities  were  computed  by  numerical  differentiation  of 
dipole  derivatives  with  respect  to  the  electric  field. 

2.3.3  Temperature  dependence  of  Raman  spectrum 

Figure  4  (a)  shows  the  Raman  spectra  collected  from  a  eutectic  Li2C03-Na2C03  (52  mol% 
Li2C03-Na2C03)  melt  over  a  band  range  of  650-1,850  cm'1  as  a  function  of  temperatures  in  a 
pure  C02  atmosphere.  For  the  carbonate  in  solid  state,  the  Raman  spectra  are  seen  to  contain  four 
basic  vibrational  modes  relevant  to  CO32'  ions.  The  two  bands  at  1,078  and  1,094  cm'1  are 
assigned  to  Vi  of  symmetric  stretching  vibrations  in  Li2C03  and  Na2CC>3,  whereas  the  bands  at 
792  cm'1  and  870  cm"1  are  assigned  to  v2  of  out-of-plane  bending  vibrations.  For  the  isolated 
CO32'  that  has  a  D3h  symmetry,  v2  is  Raman  inactive.  However,  it  is  likely  that  v2  becomes 
Raman  active  for  (Li/Na)2C03  due  to  the  distortion  of  the  CO32'  structure  imposed  by  the 
cations.  The  observed  bands  at  704  and  728  cm'1  are  assigned  to  V4  of  in-plane  bending 
vibrations.  This  mode  is  a  double  degeneration  for  the  distorted  CO32'  induced  by  Li+  and  Na+. 
As  the  CO32'  group  becomes  distorted  from  its  regular  planer  symmetry,  this  mode  splits  into  two 
components.  The  bands  at  1,375  cm"1,  1,404  cm"1,  1,531  cm'1  and  1,563  cm'1  are  attributed  to  a 
split  V3  of  the  asymmetric  stretching  vibrations  caused  by  the  existence  of  Li+  and  Na+  around  the 
CO32"  ions. 
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Figure  4  (a)  The  Raman  spectra  of 
(Li/Na)2CC>3  in  the  region  of  650-1850 
cm'1  as  a  function  of  temperature  in 
CO2  atmosphere;  (b)  Magnified  view 
in  the  region  of  1,000-1,150  cm'1 


The  shift  in  Raman  band  with  temperature  is  better  viewed  in  Figure  4  (b),  a  magnified  spectrum 
showing  the  region  of  the  major  vi-bands  at  1,078  cm'1  and  1,094  cm'1.  As  the  temperature 
increases,  the  vi-bands  for  Li2CC>3  and  Na2CC>3  shift  toward  lower  wavenumber  and  eventually 
merge  into  one  broad  peak  at  the  melting  temperature  of  490°C.  This  shift  is  a  direct  result  of 
lowered  force-constant,  elongated  C-O  bond  length  and  weakened  Li(Na)-C-0  bond  strength  by 
increasing  temperature.  When  the  temperature  reaches  the  melting  point,  the  overtone  of  the  out- 
of-plane  bending  mode  (2*V2)  appears  at  1,762  cm'1.  It  should  be  noted  that  a  broad  but  small 
peak  near  970  cm"1  only  observed  in  solid-state  carbonates  appears  not  directly  related  to  CO32 
identification  of  which  is  not  possible  for  this  study.  We  speculate  that  the  impurities  in  the 
sample  could  be  a  source  for  this  unknown  peak.  A  detailed  assignment  of  the  measured  Raman 
bands  at  different  temperatures  and  atmospheres  is  summarized  in  Table  1. 

The  most  distinguishable  features  of  Raman  spectra  in  Figure  4  are  observed  when  the  carbonate 
is  in  a  molten  state.  The  four  bands  corresponding  to  the  V3  mode  disappear  from  the  spectrum 
while  two  new  broad  bands  at  1,317  cm'1  and  1,582  cm'1  emerge  within  the  same  band  width.  A 
natural  question  to  ask  is:  are  these  newly  emerged  Raman  bands  associated  with  the  C2O52" 
species? 

Table  1  The  Raman  frequencies  measured  and  assigned  for  the  eutectic  (Li/Na)2CC>3  at 
selected  temperatures  and  atmospheres 


Atmospheres 

C02 

Air 

n2 

Temperatures  (°C) 

RT 

350 

455 

490  (melted) 

490  (melted) 

490  (melted) 

Frequencies  (cm1) 

704 

706 

706 

707 

707 

707 

1)4  (weak) 

728 

720 

715 

- 

- 

- 

792 

790 

789 

790 

790 

790 

1)2  (weak) 

879 

879 

881 

882 

885 

885 

1,078 

1,077 

1,077 

1,070 

1,072 

1,072 

Di  (strong) 

1,094 

1,090 

1,087 

- 

- 

- 
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1,375 

1,385 

1,386 

- 

1,391 

1,392 

1,404 

1,415 

1,416 

- 

1,421 

1,421 

D3  (weak) 

1,531 

1,520 

1,520 

- 

- 

- 

1,563 

1,555 

1,545 

- 

- 

- 

New  peaks  (broad) 

- 

- 

- 

1,317 

- 

- 

1,582 

Overtone(2*D2) 

- 

- 

- 

1,762 

1,762 

1,762 

Unknown 

- 

- 

972 

- 

- 

2.3.4  Atmosphere-dependence  of  Raman  spectrum 


To  answer  this  question,  we  first  measured  Raman  spectra  in  different  atmospheres.  According 
to  the  enabling  electrochemical  reaction  shown  in  reaction  (3),  the  formation  of  C2O52'  requires  a 
source  of  CO2.  Figure  5  compares  the  Raman  spectra  measured  in  N2,  Air  and  pure  CO2 
atmospheres  at  490°C  where  the  carbonate  is  in  a  molten  state.  The  bands  at  1,072  cm'1,  790  cm'1 
and  885  cm'1,  707  cm"1,  1,391  cm'1  and  1,421  cm'1,  1,762  cm'1  shown  in  Figure  5  (a)  recorded  in 
N2  atmosphere  correspond  to  the  symmetric  stretching  ( v  1 ) ,  out-of-plane  bending  (V2),  in-plane 
bending  (v4),  asymmetric  stretching  (v3),  and  the  overtone  of  the  out-of-plane  bending  mode 
(2*V2)  vibrations,  respectively.  The  Raman  spectrum  in  air  is  almost  identical  to  that  collected  in 
N2.  However,  the  peaks  at  1,317  cm'1  and  1,582  cm"1  shown  in  Fig.14  (c)  are  only  observable  in 
the  CO2  atmosphere. 


Figure  5:  Raman  Spectra  of  CO2  in  Different  Atmosphere 
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The  strong  C02-dependence  of  the  bands  at  1,317  cm'1  and  1,582  cm'1  provides  a  crucial  hint 
for  the  formation  of  C2O52'  via  the  CO2  chemisorption  reaction  C02+C032'=C2052'. 

2.3.5  DFT  Modeling 

The  theoretical  support  to  the  formation  of  C2O52"  is  provided  by  the  DFT  calculations.  The 
calculations  indicate  that  high  Raman  activities  of  C2O52'  are  within  a  band  width  of  1,200-1,600 
cm'1.  Specifically,  the  active  Raman  bands  are  predicted  at  1,366  cm"1,  1,531  cm"1  and  1,566  cm" 
1  for  Li2C205  and  1,345  cm'1,  1,547  cm'1  and  1,579  cm'1  for  Na2C20s,  respectively.  The 
recorded  spectrum  in  the  band  width  of  1,200-1,650  cm'1  from  the  MC  within  490-525°C  in  CO2 
atmosphere  could,  therefore,  be  an  overlap  of  these  characteristic  Raman  peaks  of  Li2C205  and 
Na2C2C>5  in  this  band  region.  To  deconvolute  the  two  unique  broad  peaks  around  1,317  cm"1  and 
1,582  cm'1,  we  used  Guassian-Lorentizian  function  with  the  six  theoretical  Raman  frequencies  as 
the  standards.  The  results  are  shown  in  Figure  6,  where  the  black  and  red  lines  represent  the 
measured  and  modeled  spectra,  respectively.  Also  shown  are  the  individual  spectrum  calculated 
for  the  pure  Li2C2C>5  and  Na2C2C>5,  represented  by  pink  and  blue  lines,  respectively.  It  appears 
that  the  modeled  spectrum  is  dominated  by  the  Na2COs;  only  one  peak  at  1,566  cm'1  is  visible 
for  the  O2CO3  while  the  other  two  are  too  weak  to  be  seen.  This  is  primarily  due  to  the 
differences  in  size  and  polarizability  of  the  Na+  and  Li+  ions.  Overall,  the  measured  Raman 
spectrum  is  in  excellent  agreement  with  the  DFT  calculations  if  the  temperature  effect  is  factored 
in. 


Figure  6:  Deconvolution  of  the  two  broad  bands  observed  at  1,317  cm'1  and  1,582  cm'1;  Inset:  the 
measured  and  DFT-modeled  Raman  spectra  in  a  band  width  of  1,200  -  1,650  cm'1 

In  summary,  we  demonstrate  the  first  experimental  evidence  for  the  existence  of  pyrocarbonate 
C2O52'  species  in  a  eutectic  Li2C03-Na2C03  melt  exposed  to  CO2  atmosphere  through  a 
combined  “DFT  modeling”  and  “Raman  Spectroscopy”  approach.  The  existence  of  C2O52' 
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species  is  a  key  support  to  a  new  bi-ionic  transport  model  established  to  elucidate  the 
fundamentals  of  the  high-flux  CO2  transport  phenomenon  observed  in  the  superior  mixed  oxide- 
ion  and  carbonate-ion  conducting  CO2  separation  membranes.  The  broad  Raman  bands  centered 
at  1,31V  cm'1  and  1,582  cm'1  are  a  characteristic  of  C2O52'  species  in  molten  carbonates  exposed 
CO2  atmosphere.  The  measured  characteristic  Raman  frequencies  of  C2O52'  are  in  an  excellent 
agreement  with  the  DFT-model  consisting  of  six  overlapping  individual  theoretical  bands 
calculated  from  Li2C205  and  Na2C20s. 


3.  Future  Work 

The  work  from  the  current  period  has  concluded  our  DFT  work  on  the  oxygen  reduction  in  MC 
and  the  research  will  move  to  next  stage  using  the  configuration  interaction  embedding  method. 
The  DFT  methods  have  uncertainties  when  a  charged  systems  is  considered,  evidenced  by  the 
test  case  of  oxygen  on  Ag(100).  Further  studies  will  be  demanded  to  examine  the  dependency  of 
DFT  hybrid  functionals.  Even  there  are  many  functionals  available  to  choose,  but  lack  of  guide 
for  choosing  them  accordingly  to  a  specific  system. 

For  oxygen  diffusion  in  MC,  the  DFT  study  is  the  first  step  and  provides  obtained  geometries  for 
the  Cl  studies,  which  is  hard  to  achieve  in  Cl.  We  will  use  the  optimized  geometry  from  DFT  in 
the  Cl  calculations. 
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